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ABSTRACT 


A small, uncooled liquid rocket motor was utilized to 
examine the effects of the magnitude and direction of a 
steady vortex velocity on combustion instability. The rocket 
operated on normal heptane and air, with air flow rates con- 
trollable to allow for various amounts of tangential momentum. 

Results of this investigation indicated that an energy 
feedback mechanism, through tangential momentum, did exist 
for both tangential and radial modes of instability. It was 
further shown that limit values of tangential velocity did 
exist, for a particular combustor, above and below which 


tiastability did not occur. 
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lL. INTRODUCTION 


Liquid propellant rocket engines almost inevitably 
encounter the problem of combustion instability at some time 
during engine development. In the past two decades a number 
of combustion instability models have emerged in an attempt to 
accurately predict instability boundaries. 

Combustion instability refers to pressure oscillations 
in the combustion chamber which, when encountered in a rocket 
engine, can cause an increaSe in heat transfer rates up to 
ten-fold resulting in engine destruction in the order of 
milliseconds. 

There are three principal types of combustion instability 
{Refs. l and 2]: low, intermediate and high-frequency. Low 
frequency instability, or chugging, generally refers to oscil- 
lations up to approximately 100 hertz and is the result of 
coupling between the propellant feed system and the combustion 
process. Intermediate frequency instability, typically from 
100 to 1000 hertz, is not fully understood but is thought to 
involve structural vibrations, feed system line resonances, 
vortex formation around sharp corners and high heat transfer 
rates. Pressure variations in this frequency range of insta- 
bility, reach about five percent of the steady state chamber 
pressure and tend to be more of a nuisance than damaging to 
the engine. High frequency instabilities are associated with 
the acoustic modes of the combustion gases and chamber and 


are the result of disturbances which are driven by interaction 
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of the combustion process and the chamber geometry. Of these 
instabilities, the high frequency type is the most destruc- 
tive, since it can result in pressure fluctuations exceeding 
100 percent of the steady state chamber pressure. 

The destructive nature of combustion instability ina 
rocket engine has resulted in a continuing effort toward its 
recognition and elimination. 

Analytical models have served satisfactorily to indicate 
qualitative trends in the combustion process, Bee to date 
none have been able to give quantitative information which 
would enable the prediction of the onset of instability ina 
particular rocket motor. In fact, no general agreement has 
been reached concerning the specific cause(s) of combustion 
ieeabi ities in liquid propellant rockets. 

Droplet vaporization has generally been considered to be 
the rate controlling mechanism in many liquid propellant 
rocket combustors and some stability models have been based 
upon this phase of the overall combustion process. However, 
Strahle [Ref. 3] has suggested that while vaporization may 
be important, it cannot be the major supporting mechanism in 
combustion instability. Other data [Ref. 4] also appears to 
Support this contention. 

Vortex velocity refers to the circumferential or tangen- 
tial velocity component within a cylindrical combustor. 
References 5, 6, 7, and other studies have shown that for a 
droplet vaporization limited system, the existence of a steady 


vortex velocity causes a dramatic change in the boundary 
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between stable and unstable combustion as shown in Fig. l. 

Reference 8 investigated the effect of tangential gas 
velocity as a means of controlling instability in a vapor- 
1zation-limited combustion process. Both experimental and 
theoretical evidence are presented which indicate that a 
Steady tangential velocity of low magnitude can readily excite 
the traveling (i.e., Spinning transverse) mode of combustion 
instability in a planar two-dimensional combustor. The effect 
of injecting gaseous nitrogen into a two-dimensional liquid 
oxygen, gaseous hydrogen combustor was investigated in Refs. 
9 and 10. These studies also found that transverse combustion 
instability was induced with tangential injection of nitrogen 
in the combustion chamber. 

An initial investigation was conducted to study the effect 
of a steady vortex gas motion on combustion stability in a 
small uncooled three-dimensional research rocket motor [Ref. 
11}. The specific objective of this investigation was to ex- 
tend the work of Ref. 11 in order to determine the effects of 
Magnitude and direction of vortex flow velocity on combustion 


instability in a rocket combustor of practical configuration. 
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II. GENERAL DISCUSSION OF THE PROBLEM 


In general, Raleigh's criterion [Ref. 12] states that 
energy release in phase with a wave tends to drive the wave, 
and out of phase tends to dampen it. As mentioned above, 
several models, based on this criterion, have been proposed 
in the literature which can be used to predict the effect of 
a steady vortex gas flow on combustion instability. 

Priem developed a model [Ref. 5] which considered the 
various physical processes (droplet atomization, vaporization, 
chemical reaction, injection, etc.,) in determining combustion 
instability. By assuming that the combustion rate was con- 
trolled by the slowest process, he was able to theoretically 
predict a Minimum pressure disturbance that would be amplified 
by the combustion process. 

The model employed was a two-dimensional annular section 
with small thickness and length. The transport equations 
utilized in Priem's mathematical model are developed in Ref. 
13. Several plausible assumptions are made which simplify the 
model sufficiently for computer solution. 

If droplet vaporization is considered to be rate limiting, 
an additional expression (model) is needed for this "energy 
SOurce process." One vaporization model has been developed 
by Priem in Ref. 14. Models for other “rate limiting" steps 
are presented in Ref. 5. Itis shown in Ref. 5 that vaporiza- 
tion and atomization are strongly dependent on tangential gas 


velocity. If injection or chemical kinetics were rate limiting, 
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the energy source process is shown to be independent of 
tangential gas velocity. 

If it is assumed that vaporization is rate controlling, 
Priem's model, with no vortex flow, predicts stability bound- 
aries as depicted in Fig. 1 from Ref. 5. To determine the 
Significance of vortex flow in stability, calculations were 
performed using the vaporization model with a superimposed 
constant vortex velocity. Results are ou inshag. {irom 
Ref. 5) where it is seen that even a small vortex velocity 
may effect the stability boundaries dramatically. Although 
experimental evidence supports the qualitative trends shown 
in Fig. 1 and 2, disturbance levels for onset of instability 
are found to be one to two orders of magnitude higher than 
predicted [Refs. 5, 6, 8, and 9]. 

Heidmann [Refs. 8 and 9) examined in some detail the ef- 
fects of vortex flow in a two-dimensional circular liquid 
propellant combustor. These studies involved the injection of 
a tangential flow of an inert gas into a stable combustor. 
variable amounts of nitrogen were introduced in order to find 
the magnitude of Vo required to produce instability. Heidmann 
found that the combustor, which was stable under conditions of 
no tangential gas flow, had definite stability boundaries 
dependent on magnitude of tangential gas velocity. He also 
found that combustor pressure oscillation amplitude was nearly 
linearly related to tangential gas velocity. 

Heidman and Feiler [Ref. 8] proposed a model to explain 
the experimental and theoretical evidence that a tangential 
velocity did excite unstable burning. In this analytical 
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model it was assumed that the energy released during combus- 
tion depended on the vaporization rate which varied as 


(p(t) [av(t) |) 2/2 


where [AV(t)| represents the net velocity dif- 
ference between the liquid and gas. The extent of driving or 
damping the instability was determined by comparing the time 
variation of energy release with pressure variation and com- 
puting the net in phase or out of phase energy addition. This 
can be more readily visualized if the energy release is con- 
sidered to vary as |v|-rather than (olav]) 7/2. In this case 
the energy/pressure relationships can be explained in a 
Simplified manner using Fig. 3 from Ref. 8. Figure 3(a) de- 
picts the pressure variation with time for a traveling trans- 
verse wave with no steady tangential velocity component. 

Figure 3(b) shows a similar velocity variation. If the energy 
release is linearly related to the absolute velocity, the 
energy release can be depicted as in Fig. 3(b). It is evident 
that energy release is symmetric reaching a peak twice per cycle. 
Therefore, according to Raleigh's criteria, no energy is added 
tending either to drive or dampen the wave. If now a steady 
velocity in the plane of the wave is superimposed as shown in 
Fig. 3(c), the in-phase energy release increases at the expense 
of the out of phase energy release. Thus, the wave 1s driven. 
If the velocity had been biased in a direction opposite to the 
Wave motion, the net energy addition would have been out of 


phase with the wave,, and the wave would have been dampened 


(Fig. 3(d)). 
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Nearly all the theoretical and experimental work on the 
effect of tangential gas velocity on combustion instability 
to this time has been two dimensional, and it was the purpose 
of this investigation to address the realistic problem of a 
three-dimensional combustor. 

Kiels initial investigation [Ref. 11] attempted to show 
whether the evidence found for a two-dimensional combustor 
were also applicable to a three-dimensional rocket. He found 
that with a first tangential instability [Ref. 2] a moderate 
swirl in one direction tended to decrease the instability 
intensity whereas flow in the opposite direction tended to in- 
crease it. For small amounts of swirl, there was virtually 
no effect on instability peor Gt ee) eo the direction of 
Swirl was changed, as was also the case for a high degree of 
Swirl. However, instrumentation precluded determination of 
instability direction or accurate tangential velocity 


measurement. 
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EIT. EXPEREMENTALE APE ROACH 


A small uncooled rocket motor, described in Section IV, 
was used in conducting this investigation. Normal heptane 
and air were burned in the combustion chamber with mixture 
ratios from 12.1 to 15.2 at a nominal chamber pressure of 
260 psia. 

Air flow was divided into primary and secondary flow 
regimes in which secondary flow direction could be varied from 
the control console in a continuous test from clockwise to 
counterclockwise to straight. The manner of accomplishing 
this will be described later. 

Secondary air mass flow rates were varied from 10 to 25 
percent of the total air mass flow rate. 

The motor was instrumentated with high frequency response 
pressure transducers in order to experimentally determine the 
effects of both the magnitude and direction of tangential gas 
velocity on the magnitude and direction of the combustion in- 
Stability. Insert elements (Fig. 4), which directed secondary 
air flow, were varied in order to vary the magnitude of Vg on 
different tests while maintaining fixed secondary mass flow 
mace . 

A typical firing sequence will be discussed below. First, 
power was supplied to the igniter whereupon oxygen and methane 
were injected, mixed and ignited. Then secondary air was 
allowed to flow first in a clockwise or counterclockwise 


direction until a steady state flow rate was achieved. Primary 
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air flow was then initiated and after a few seconds fuel was 
injected at which point ignition took place. This configura- 
tion was allowed to burn for one to two seconds whereupon 
Opposite swirl was selected followed by straight secondary 
air flow. Each flow regime was allowed to burn for one to 
two seconds. After a four to six second total burn, the 
fuel was cut off, followed by primary air, allowing the sec- 
ondary air to continue to purge and cool the motor. Visual 
inspection of the motor was conducted before and after each 
run to ensure system components were intact and operational. 
Prior to each day's run instrumentation was allowed to warm 
up for two hours after which strain gage transducers were 
zeroed. All pressure transducers were calibrated (over full 
utilized span) with an Ampthor dead weight tester. 

The magnitude and direction of the secondary air velocity 
could be estimated at the injector face from measured temper- 
atures and mass flow rates. However, in the stability models 
it is not the inlet V, that is of importance, but rather the 


0 


Magnitude of V, of the total gas flow within the combustor. 


9 
Determination of tangential velocity in the combustor during a 
hot firing presents some difficulties. A simplified approach 
was used in this study. Chamber gas velocity magnitude and 
direction were determined in cold flow tests and then scaled 
to the high temperature and pressure of actual test conditions 
in a method described in Section V. Although this procedure 


may not provide exact data for the hot firing condition, it 


. does provide a means for accurately comparing the effects 
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of various magnitudes and directions of secondary air flow 
on the magnitude and direction of the velocity of the total 
gas flow. 

Gas velocity was determined as described below. Air alone 
was allowed to flow within the combustor in both the primary 
and Secondary regimes. A three hole pressure probe (shown in 
Fig. 5), was used to determine the flow direction at various 
radial positions. Once known, a velocity profile was deter- 
Mined, and a Swirl component ascertained. These tests were 
conducted for various percentages of secondary air flow. 

In order to predict the instability frequencies which 
‘Might be encountered in the rocket engine, it is necessary to 
solve the wave equation which governs the flow field in an 
acoustic medium [Ref. 15]. Although the apriori application of 
the wave equation is questionable since a small disturbance 
assumption is made, it has been found empirically that the 
solution 1S valid even for oscillations of large amplitude 
[Ref. 16]. 

Solution of the wave equation yields frequencies that are 
a function of combustor radius, length, speed of sound in the 
Medium and roots of a Bessel function of order N [Ref. 15]. 
For the combustor used in this investigation, the frequencies 


were determined to be as follows: 


Instability Frequency in Hertz 


Tangential Radial Longitudinal 
lst mode 5280 i, ood 3000 
2nd mode 8750 21,000 6000 


rg 








‘Identification of the type of instability could be deter- 


Mined not only from frequency data but also phase relationships 


and characteristic wave form. Mixed modes of instability can 


also exist but the predicted frequencies are far above any 


encountered in this investigation. 
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IV. DESCRIPTION OF APPARATUS 


A. GENERAL 

A small uncooled research rocket motor was constructed in 
the main by Kiel [Ref. 11] for the purpose of studying com- 
bustion instability in a biphase (liquid-gas) rocket combus- 
tion chamber. Both primary and secondary air flow metering 
devices were strengthened and redesigned to more accurately 
control and measure air flow rates. A new fuel injector (to 
be described later) was designed to provide a more uniform 
distribution and better atomization of fuel. It was felt that 
this would improve combustion efficiency (n.,) in all oper- 
ating regimes of burning. High frequency response pressure 
transducers were added in the combustion chamber to allow 
determination of instability type and direction of travel. 

The rocket motor was mounted on a thrust stand located in 
the rocket test cell (Fig. 6). Rocket operation was directed 
from the control room from two control panels (Fig. 7a and b) 
overlooking the test cell. All data was recorded on instru- 


mentation location in the control room (Fig. 8). 


Bee ROCKET MOTOR 

Ehe tocket motor isi shewn in Fig. 9. ft is cylindrical 
in shape and made in several sections in order to allow for 
expansion and modifications. Constructed of stainless steel, 
the rocket motor had a contraction ratio of 9.4. The nozzle 


(Fig. 10) was bolted separately to the combustion chamber to 
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allow for future modifications in contraction ratio or nozzle 
geometry. Two static pressure taps were located in the plane 
of the igniter and just upstream of the nozzle. High frequency 
response pressure transducers (Fig. 11) were located at 90- 
degree intervals half way down the combustion chamber. Two 
Kistler (Model 615) and a photocon 352A pressure transducer 
were used in order to determine instability frequency, inten- 
Sity and direction of travel. A photo port (Fig. 11) was 
available in the same plane as the high frequency pressure 


transducers for streak photography observation. 


eeekUil SUPPIEY SYSTEM 

A schematic of the fuel supply system is shown in Fig. 12 
and photographed in Fig. 13. Fuel was supplied from a 20- 
gallon tank located in a remote cell. The fuel tank could be 
pressurized with nitrogen to 1000 psi. A liquid level indi- 
cator on the fuel tank, and one on the fuel control console 
(Fig. 14), indicated fuel quantity in percent of tank capacity. 
Fuel passed through a filter and an emergency shutoff valve 
activated from control console A (Fig. 7a). It then passed 
through a microporous filter and a Jamesbury on-off actuator 
eentrolled from console B (Fig. 7b). Locally manufactured 
Ccavitating venturis were employed to control the fuel flow 
rate. The venturis were calibrated by measuring flow rate as 
a function of upstream pressure. Upper limits of back pres- 
Sure were established at which point cavitation ceased to 
occur. Fuel manifold pressures exceeding these limits would 


allow undesirable interaction of the feed system with the 
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combustion process and were to be avoided. Fuel then passed 
into the primary fuel injector manifold (Fig. 15) and then in 
order to pass through the labyrinth of primary and secondary 
Eipiag Sep ananmuas routed to a secondary fuel manifold just up- 
stream of the injector face. The centrally located fuel 
injector tube had sixteen circumferentially located holes 
drilled back toward the injector face at a 45-degree angle to 
provide for better atomization and mixing. In addition, four 
injectors elements were located at 90-degree intervals around 
the center injector element on a 3-1/2 inch diameter circle 


with single centered holes to ensure complete fuel distribu- 


€1i0n within the combustion chamber. 


D. AIR SUPPLY SYSTEM 

The air supply system is shown schematically in Fig. 12. 
Air was supplied by two banks of air tanks which were capable 
of supplying air pressure up to 3000 psi. The tanks were 
pressurized by an Ingersoll Rand 3500 psig capacity air com- 
pressor. Air pressure to the motor was regulated through a 
Grove powreactor dome controller. Air flow was divided into 
primary and secondary regimes where flow rates were measured 
using thin plate, sharp-edged orifices with 1/2d and ld pres- 
Sure taps, installed according to ASME specifications [Ref. 17]. 
Inlet air temperature was measured with a copper-constantan 
thermocouple located downstream of the orifice plates. Primary 
air was injected into the combustion chamber around the cen- 
trally located fuel injector [Fig. 18] and circumferentially 


located injector elements. 
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The center primary air inlet to the combustion chamber 
was configured such that the air flow was choked at the air 
manifold exit and then expanded to supersonic flow before 
entering the combustion chamber. Upon entering the combustion 
chamber, the flow shocked cauSing mixing of fuel and air at 
the injector face. The primary purpose of the airflow around 
the circumferentially located injectors was fuel atomization 
and its flow rate was dictated by the desired droplet size 
[Ref. 4]. Secondary air was directed into the secondary air 
manifold (shown schematically in Fig. 16 and pictorially in 
Fig. 17) where it was possible to select any combination of 
the three flow regimes: clockwise, counterclockwise, or 
Straight. Each airflow regime was injected through four tubes 
Spaced at 90-degree intervals around the injector face (Fig. 
4). The straight-through mode passed directly into the com- 
bustor. The clockwise and counterclockwise flows were de- 
flected at angles of 20, 30, or 40 degrees at the injector face 
by brass inserts held in place between the air orifice and air 
swirling plates (Fig. 4). The lower mass flow rate of second- 
ary air was restricted by the ability of the flow orifice to 
properly meter the flow, and the upper limit by the amount of 
flow the secondary air solenoid valves would pass. It was 
found that low air mass flow rates produced an oscillation 
within the air feed system that could not be eliminated, thus 
restricting the allowable flow rates. Numerous attempts were 
made to isolate the problem, as it was not known initially 
whether the oscillations were due to instrumentation or due 


to the propellant supply system. Various air flow metering 


24 








orifices were used within the limitations prescribed by the 
ASME power test code [Ref. 171, But ta neowavalierina wan 
secondary flow regimes were run individually with all instru- 
mentation powered by separate power supplies with flow both 
directed to the rocket motor and fed directly to the atmos- 
phere. This showed that the problem was not with the instru- 
mentation or the downstream configuration. It was then 
concluded that the resonance was within the air feed system 


but efforts failed to eliminate it. 


me LGNITION 

The ignition system is shown schematically in Fig. 19. 
The combustion mixture was ignited with a methane oxygen torch 
shown schematically in Fig. 20 and photographed in Fig. 21. 
Current was supplied to a matrix of four nichrome wires 
insulated by a mullite tube. Gaseous methane and oxygen were 
injected separately in concentric tubes around the mullite 
and mixed for a short distance prior to exposure to the spark. 
tiie igniter continued to burn until de-activated by a pressure 
sensing switch which shut off the methane, oxygen, and spark 
when the chamber pressure reached about 50 percent of its ex- 


pected steady state value (Fig. 22). 


FF. INSTRUMENTATION 


A summary of the instrumentation used is listed in Table 


Primary and secondary air pressures and differential pres- 
Sures were used with temperature readings to determine air 


flow rates. Air manifold pressure was required to determine 
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possible instability coupling) wath the alr tecau>s oon 
Steady air manifold pressure was necessary to ensure no 
coupling. Venturi pressure allowed for the determination of 
fuel flow rate and it was necessary that fuel manifold pres- 
sure be below a predetermined value to ensure that cavitation 
occurred in the venturi. Pressure readings from the low 
output strain gage transducers, were amplified by a factor of 
100 and then fed into the control room where the signals were 
attenuated (Fig. 23) and recorded on a Honeywell 1508 Visi- 
corder oscillo-graph (Fig. 25). Higher output signals were 
attenuated and recorded directly on the visicorder. High 
frequency response pressure transducers were utilized for in- 
Stability detection within the combustion chamber. These 
transducers were located in the combustion chamber (Fig. 5) so 
as to give maximum information about the instability (direc- 
tion of travel, type, strength, and frequency). The pressure 
variations from the photocon and Kistler transducers were 
recorded as an Ampex CP-100 magnetic tape recorder at high 
tape speed. The signals on the Ampex were then played back 
at slow speed onto the visicorder operating at its fastest 
speed (80-IPS). This procedure was necessary in order to view 
the three traces for the entire run simultaneously where fre- 
quency, phase and amplitude could be determined. The three 
Signals were also photographed simultaneously on a three 
channel oscilloscope for phase relationship determination. 
The purpose for examining these traces on the oscilloscope 


was to ensure that no phase information was lost in 
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transferring the high frequency output from the tape to the 
visicorder. Examining these traces on the oscilloscope was 
accomplished on the first few runs only until it had been 


shown that no phase information had been lost. 
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VV. JRESULES 


A summary of the data obtained from this investigation 
is presented in Table II. Only data from successful firings 
where instrumentation and rocket motor were operating satis- 
factorily are reported. Run Numbers 1 through 6 employed 
20-degree inserts and runs 7 through 17 employed 40-degree 
inserts. 

Fuel flow, ane flow rates and mixture ratio were computed 
using a simple data reduction program. Average chamber pres- 
sure between the two chamber pressure transducers is reported. 
It should be noted that disparity was never greater than 
2 percent between these two transducers. 

Difficulty was encountered in obtaining rocket ignition 
during tests at high secondary air flow rates with the 40- 
degree inserts. Altering the secondary air flow sequencing 
resulted in successful ignition. For this reason, sequencing 
was altered in almost all runs utilizing the 40-degree inserts. 
c* efficiency was calculated for each run using a theoretical 
c* predicted by using the computer program of Ref. 18. 

Theoretical c* values are shown in Fig. 24. Instability 
frequency and direction (when determined) are also presented 
in Table II. In addition, calculated values of the tangen- 
tial velocity at the wall and the ratio of injected tangential 


to injected total momentum are included. 
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One of the parameters of great importance in this anal- 
ysis, was tangential gas velocity. Ideally it would be desir- 
able to obtain this information from an actual firing. Due to 
the excessive gas temperatures within the combustion chamber, 
this was impossible. It was decided to obtain this information 
by scaling data obtained in a cold flow situation, with mass 
flow rates equivalent to those used in hot firings. This was 
accomplished as follows: A 3-hole pitot probe was rigged in 
the same plane as the high frequency pressure instrumentation. 
The radial position of the probe was varied in one-fourth inch 
increments across the combustion chamber to determine a flow 
direction and velocity profile. This data was obtained at 
0, 10, and 22 percent secondary mass flow rates, with both the 
20- and 40-degree inserts. Velocity values for secondary flow 
rate between these values were obtained by linear interpolation. 
In order to estimate the tangential gas velocity in the hot 
firings, the following analysis was employed. The mass flow 
rate during a hot firing was almost equivalent to that during 


a cold firing (neglecting the small mass flow rate of fuel) 


(ODN ot oO timend 1) 
J m8 Oley 
(apAV) not ~ (RrAY) cola 


Pressure (P nola? : velocity Oe al i and temperature (T ola? 


were determined for each cold flow test. For actual rocket 
firings, pressure (PL ot) was measured and temperature (Ty ot? 


was calculated using the experimental values of mixture ratio 
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and combustion efficiency. Assuming that Boat = A ola! 


Eq. la could then be used to calculate Vise a Assuming that 
the velocity direction did not change between hot and cold 
flow tests, Vo for the hot firings could be calculated. It 
was also assumed that the magnitude of any steady tangential 
gas velocity induced by the instability itself, was negligible 
with respect to the injection induced Vo- It 1S recognized 
that although this analysis is not precise, it should give a 
good estimate of the velocity values in a hot run. 

While obtaining the cold flow data, it was noticed that 
there existed regions of reverse flow close to the walls. 
Since it was not possible to obtain reliable velocity informa- 
tion very near the walls, it was decided to verify this using 
a tuft analysis. This analysis confirmed the reverse flow 
regimes and large amounts of recirculation in the chamber. It 
was hoped that the pitot probe results could be substantiated 
through conservation of momentum by comparing inlet momentum 
at the injector face with measured momentum at the plane of 
the pitot tube. The reverse flow near the wall precluded this 
analysis. 

The momentum ratio of tangential momentum to the total 
momentum was computed at the injector face through knowledge of 
the primary and secondary flow rates, direction, temperature 
and pressure. 

The results of the cold flow analysis are shown in Figs. 
25 and 26. With the velocity magnitude and direction known, 


it is possible to determine the components of velocity in the 
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axial and tangential direction. These components are non- 


dimensionalized by ‘o which occurred along the chamber 


ax 
centerline. 

A typical visicorder output is shown in Fig. 27. This 
pressure-time trace was used to determine various output 
pressures and finally determine flow rates, and engine per- 
formance parameters. High frequency data output from the 
Kistler and Photocon transducers are shown in Fig. 28. The 
Same output is shown in Fig. 29 with an expanded time scale. 
The output in Fig. 28 was used to determine the instability 


persistance, while the output in Fig. 29 was used to determine 


phase and amplitude information. 
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VE. DESCUSSTON OF RESULTS 


The instabilities encountered in this investigation were 
first radial, first tangential and second tangential. Type 
of instability was identified through phase and amplitude 
relationships obtained from the three transducer traces. 
Most of the instabilities encountered were first radial and 
thus not helpful in demonstrating the original objectives. 
Second tangential instabilities were more useful but deter- 
mination of direction of travel was not possible. Spinning 
first tangential instabilities where direction of travel was 
readily identifiable would have been most conducive to study- 
ing the objectives of this study. These occurred only in 
runs 9 and 13. Run 13 was spontaneously unstable with a first 
radial instability during counterclockwise swirl operation. 
Upon selecting clockwise swirl, the instability switched to 
a spinning first tangential traveling in the opposite direction 
to the swirl. Instability amplitude during the counterclock- 
wise phase was about 33 percent of that in the clockwise phase, 
which is consistent with Heidmann and Feilers hypothesis that 
energy added out of phase with a traveling wave tends to 
. attenuate it and in phase tends to drive it. Attempts to 
repeat this test condition were unsuccessful. 

The fuel injector used by Kiel and in his study was a 
Single element injector which resulted in a fairly concentrated 


energy distribution within the combustion chamber. With this 
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fuel injector it was not possible to support combustion in 

the straight through mode of secondary air flow. A multi- 
element fuel injector (Fig. 18) was incorporated in the 

rocket motor for this investigation in an attempt to estab- 
lish efficient combustion during all modes of operation. 

That is, 1t was desired to ensure equally complete combustion 
during burning with and without vortex flow. It was found 
that the straight through air flow mode of combustion could be 
sustained with the multi-element injector. However, the vor- 
tex flow continued to cauSe more efficient burning due to 
better droplet atomization and increased fuel residence time 
within the combustor. This resulted in combustion efficien- 
emes with vortex flow of from .95 to 1.10 and from .80 to .87 
with no vortex flow. This causes a slight degradation in the 
results since the effect of combustion efficiency and tangen- 
tial gas velocity cannot be uncoupled in the results. There 
were two notable exceptions to the combustion efficiency 
norms. On run 6, a combustion efficiency of 1.01 was attained 
in the "no swirl" mode. This resulted from loss of the center 
fuel injector during mid-run. Loss of the center fuel element 
caused fuel and air to be premixed in the air manifold with 
subsequent increased combustion BEeiCHenere Fuel injector 
loss was thought to have occurred during the straight through 
flow because a drop in air manifold pressure was observed at 
that time. Higher than usual ae in the straight mode also 


occurred on run ll, and its occurence remains unexplained. 
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Although a calculated n., above 1.00 is not theoretical- 
ly possible, it is not unusual and is thought to be caused 
by an effective decrease in throat area due to boundary layer 
buildup and/or an effective throat reduction due to vortex 
flow. 

Nearly all runs exhibited the tendency toward higher 
amplitudes of oscillatory pressure with more efficient com- 
bustion ier (see Fig. 30). Mixture ratio did not have a 
Significant effect on the instability amplitude (see Fig. 31). 

In general, the most unstable conditions occurred in the 
clockwise mode of operation as was found by Kiel who used the 
Same basic combustor, but with a different injector. 

In those tests where a first radial instability persisted 
in both the clockwise and counterclockwise mode, the insta- 
bility amplitude was nearly constant. During the straight 
mode, the instability would disappear or switch stability 
modes to a first tangential instability with greatly reduced 
magnitude. This indicated that an energy coupling mechanism 
may exist due to tangential gas velocity. The tangential gas 
velocity appeared to drive the instability from a small ampli- 
tude first tangential mode to a large amplitude first radial 
mode. 

In run 6 the injector failure during the straight mode 
allowed fuel and air to be premixed in the air manifold. The 
first radial instability which existed in the counterclockwise 
and clockwise modes of air flow persisted in the straight mode 


but with reduced amplitude. 
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In addition, inspection of Fig. 32 reveals that the 


magnitude of V, affected the mode of instability. In general, 


Q 
tangential velocities between 45 and 55 feet per second and 
75 and 85 feet per second induced radial instabilities. 
Values of Ve above, below and between these ranges induced 
first tangential or second tangential instabilities. This re- 
sult is not predicted in the two-dimensional models. 

On several runs (not reported) in which very high tangen- 
tial velocities were present, it was not possible to ignite 
the motor due to excessive radial acceleration of the fuel- 
alr mixture. Inspection of the rocket motor after these un- 
successful firings showed fuel stream traces almost 
circumferentially around the combustion chamber. 

These results indicated that one of the primary effects 
of the tangential gas velocity was that it both increased the 
atomization and vaporization rate, which changes the axial 
distribution of available energy to a more concentrated 
profile and changes the radial distribution of available energy. 


At very high values of V. most of the fuel was concentrated near 


Q 
the wall which was conducive to tangential instabilities of 
stable operation if the motor could initially be ignited with 
this concentrated fuel distribution. 

In addition to the effects of V, on the radial mode of 


G 


instability, the magnitude and direction of V, may also affect 


Q 
tangential modes of instability. Most of the combustion 
instability models deal with tangential instabilities since 


radial instabilities are not permitted in a two-dimensional 
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model. Subsequent discussion will deal with the first tangen- 
tial and second tangential experimental results of this 
investigation and the first tangential experimental results 
reported by Kiel [Ref. ll]. 

There were two runs where a second tangential instability 
persisted during both clockwise and counterclockwise modes of 
operation. In run 12, the instability strength remained 
constant whereas in run 3 the amplitude was reduced. It is 
felt that the amplitude reduction occurred due to energy being 
added out of phase with the instability direction of travel. 
However, the spin direction of the instability could not be 
ascertained on this run. 

The weakest tangential instabilities tended to be the first 
tangential. The largest and smallest values of tangential 
velocity produced Bs values from zero to 0.4, whereas midrange 
values of Vo produced stronger instabilities (See Fig. 32). 
Calculated tangential gas velocity near the well varied between 
28 and 92 feet per second in this investigation. This suggests 
that vi limits may exist above and below which tangential veloc- 
ity does not drive the tangential instability. This result 
was also evident in Kiel's investigation. Cold flow data from 
this investigation was used to calculate tangential gas veloc- 
ity for his study. Ve calculated in this manner ranged between 
Q2 and 73 feet per second. 

In the theoretical model employed by Feiler and Heidmann 


{Ref. 8] an upper limit of Vo is predicted. The model predicts 


that, up to a point, increases in the in-phase energy release 
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(tangential gas velocity) results in stronger instabilities. 
Upon reaching a critical value of tangential velocity, further 
increases serve only to increase the average energy added to 
the system rather than the oscillating portion. High values 
of Vo may also increase atomization and vaporization rates 
such that they may not be the rate limiting step(s) in the 
combustion process [Fig. 32]. 

It is notable that on several runs with high combustion 
efficiency (runs 14, 16, 17) and low tangential velocity, the 
combustor was stable. Although runs 16 and 17 were unstable 
initially during the transient, they recovered and remained 
stable thereafter. This stability occurred at the minimum 
tangential velocity of about 35 to 45 fps., which was almost 


precisely what Heidmann reported for the minimum V. required 


Q 
to induce a tangential instability in his tests with a two- 
dimensional combustor. The linear increase in oscillatory 
pressure amplitude with increased tangential velocity as found 
by fewdriann and Feiler [Ref. 5] in the 2-D combustor could not 
be verified in this investigation because the mode of stability 
changed from tangential to radial at intermediate values of Va: 
In order to determine the magnitude of Vig required to 
induce spontaneous tangential combustion instability in 
Priem's model, it is necessary to first calculate a value for 
the burning rate parameter —— Although the average burn- 
ing rate of the propellant per inch (m) is not known for this 


combustor, it can be approximated. Since the combustion effi- 


ciency is generally about 100 percent, it might reasonably be 
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assumed that it takes the complete length of the combustion 
chamber to vaporize all the fuel. Using a combustion chamber 
length of 6-1/2 inches and assuming a constant value of "m" 
yields an "m" of .154. Then K = -033. Using Priem's insta- 
bility model results, as shown in Fig. 1 from Ref. 5, it is 


predicted that spontaneous instability would occur for V,>20 


6 
foce This is in qualitative agreement with the xesultes lof 

this investigation, although several runs were also unstable 

in the straight mode of air flow. Both Priem's and Heidmann's 
analytical models and Heidmann's experimental results show a 
threshold value of Vo above which the combustor 1S unstable and 
below which stable operation results. In general, the results 
of this investigation (Fig. 32) agree with this. However, it 
should be emphasized that neither Heidmann's experimental 
results nor Priem's analytical model predict an upper limit 

on tangential gas velocity as discussed above. In was found 

in this study that a limit value of Vo did exist above which 
the rocket was stable. 

Instead of discussing Vo effects, it may be of equal im- 
portance to consider tangential momentum which has been shown 
to be of importance in solid propellant combustion instability 
[Ref. 19]. A plot of the injected tangential momentum vs. Bs 
is shown in Fig. 33. This figure indicates a limit on the 


sensitivity of instability strength to the ratio of tangential 


to total momentum as was also found for Vo. 
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The quadrant configuration of the fuel injector 
resulted in an energy distribution that invited a second 
tangential and/or first radial instability. The unavail- 
ability of additional high frequency response transducers 
precluded determination as to whether a second tangential 
instability was spinning or stationary. For this reason, 
it was desired to create a spinning first tangential insta- 
bility as this would lend to a more definitive analysis. 

To this end, the circumferentially located fuel injectors 
were sealed to closely duplicate a fuel injector used by 
“Kiel [Ref. ll], whose investigation produced only first 
tangential instabilities. This effort was successful in 
producing a standing first tangential instability upon rocket 
ignition. However, the instability did not persist and was 
attenuated after about .l1 sec. These runs had the smallest 
amount of tangential gas velocity SI, which may well have 
-been the threshold value above which an instability would 
have persisted as discussed above. Runs 16 and 17 represent 
those firings in which both straight and counterclockwise or 
Straight and clockwise modes were fired simultaneously. 
Rocket ignition was not possible with only one mode selected 
as was also reported by Kiel. 

Notable in this investigation is the closeness of fre- 
quency for the second tangential and first radial insta- 
bilities. This result is not predicted by existing models. 
Powell and Zinn [Ref. 20] have presented a solution to the 
nonlinear wave equation in three-dimensional space, using 


the method of Galerkin. Their analysis predicts the second 
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tangential and first radial modes of instability to oscil- 
late at twice the frequency of the first tangential mode. 
The analysis further predicts first tangential instability 
frequency to be quite close to those predicted by linear 
theory. These results compare well with the results of this 
investigation. | 

In run 10 the mode of instability switched from second 
tangential to first radial in mid run during the clockwise 
mode of operation. This may also have occurred on other runs, 
but it was impossible to view the entire run so that detection 
of this phenomenon was not assured. This phenomenon of 
coupling between various instability modes is consistent with 


the results of Zinn and Powell [Ref. 20]. 
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CONCLUSIONS 


It was found in this investigation that for a three- 
dimensional ligquid-gas combustor, tangential momentum and/or 
gas velocity had a significant effect of combustion 
instability. 

First radial instabilities were unaffected by direction of 
tangential gas velocity but occurred only at specific values 
of Vg- This suggests that the tangential gas velocity (mom- 
entum) can affect the first radial instability by changing the 
radial energy distribution within the combustor. 

When second tangential instabilities were encountered 
during operation with swirl in a particular direction, changing 
the direction of swirl caused the instability to be attenuated 
or driven. 

When first tangential instabilities existed, instability 
direction could be determined. When instability direction was 
opposite to the swirl a dramatic reduction in instability 
strength resulted. 

Upper and lower limiting values of he were found above and 
below which the combustor was stable. The lower limit was 
found to be in qualitative agreement with the experimental 
results of other investigators using different combustor con- 
figurations. In addition, this lower limit was also in quali- 
tative agreement with the analytical models of Feiler and 


Heidmann [Ref. 8] and Priem [Ref. 5]. 
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Coupling between the various modes of combustion insta- 
bility (first radial, first tangential and second tangential) 
was also observed in agreement with the work of Zinn and 


Powell [Ref. 20]. 
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Figure 9. Rocket Motor 
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